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^ I Abstract 



o 



We construct the heavy baryonic currents by using the Bethe-Salpeter wave func- 
tions in the heavy quark hmit. We discuss the one- loop renormahzation of these heavy 
baryonic currents as well as their two-point correlators up to the order l/M^. For a 
special case, we do the QCD sum rule for masses of the doublet (3/2,5/2). 



^ . Many baryons with spins larger than 3/2 have been found in the experiments. A successful 

^ ! description of these particles in theory will further confirm the quark model. However, up 
O \ to now there is few systemic methods to study such kind of particles. On the one hand, any 
^ I relativistic model for a three-body system will be very difficult to be dealt with. On the 
other hand, there are also some problems to apply the non-relativistic models to baryons, not 



only because it is still difficult to get a solution of a three-body system but also because of 
suspectable validity of the non-relativistic approximation. In order to get more knowledge of 
the excited baryons, it may be useful to start our study from the simplest system. The heavy 
^1 baryons, which contain only one heavy quark, probably are what we are looking for . In such 

>-Ch ■ systems, the spin of the heavy quark decouples with the light freedoms when the mass of the 
heavy quark goes to infinity. This is so called the heavy quark symmetry (HQS) Then 
^ ! the total spin of the heavy baryons can be constructed in the frame of the diquark picture. 
^ \ The effects of HQS breaking can be taken into account perturbatively in the expansion of 
1/Mh in the frame of the heavy quark effective theory(HQET) 0. Conversely, the property 
of the excited heavy baryons could check whether HQET works well in such a system( it is 
still amazing why the life time of Af, is so smaller than that of B). Studying excited heavy 
baryons also provides useful information for semileptonic decays of heavy baryons which 
can serve to determine Kobayashi-Maskawa matrix elements. The reliable estimate of the 
widths of Afe — »• Aj:^**^'^ -|- lu, which probably are small, could give a more tight constraint on 
Afe^ A, + Zz/ i. 

Based on HQET, a non-perturbative method is still needed to deal with the heavy baryons. 
QCD sum rule as a model-independent method is very powerful and has wide applications 
0. Our discussion in this manuscript will be restricted in its framework. In the following 
we firstly discuss the choice of the baryonic currents in section 1. Then the general formulas 
of two-point correlators are discussed in section 2. In section 3 we give an example of the 
calculation for the masses of the heavy baryon doublet (3/2,5/2). 
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1 Baryonic currents 

In the heavy quark hmit, it is possible to choose a baryonic current that only annihilates the 
heavy baryon state with the same quantum number. This unique character is also valid up 
to order 1/M/j in QCD sum rule. In th practical application, this character is good enough 
at least for b heavy baryons. 

There are many ways to construct such baryonic currents. In this paper, we would like 
to deduce them from the Bethe-Salpeter wave functions. Let us consider the heavy baryon 
wave function 

-irngv-xi f J^^V_J^^^ -ip-x -ik x\ 



{Q\Th{x,)qi{x,)q,{x,)\B) = e— - j ^^^^e-^P-^e-''-^ x{v,P,k), (1) 

where x = X2 — x^, x' = Xi — ^^^^,v is the velocity of the baryon. We have already set 
m/j 00 and suppressed the color index for convenience . For an arbitrary spin n + |, 
x(f ,p, k) has two independent structures [^]. For the parity P = (—1)", 

S{lJ.l---fJ,n + l) 

and for P = , 

X{v,p,k) = ( ^ l5lTu^'-^-^')^,,...,^_^,^,C, (3) 

where u'^'^'"'^" is a general Rarita-Schwinger tensor spinor, J2s{tii---iJ.„) makes indices Hi - ■ ■ Hn 
symmetric and is a transversely symmetric tensor function of 7_l, p± and k^. The 

7_L (so as p_L and k_[_) is defined as 

7^ = - (4) 

The first and second state of (0) degenerate with the second and first state of (0) resp. n = 
is a special case, in which the first state of (0) is a singlet. 

Then the baryonic currents can be constructed straightly. For instance, we can choose 

r^--^-{x) = r'|--^"'^--^"/.(x)j,,..,„(x) (5) 

for the first state of (^, where 

ln+ I s{l-n) 

and where 



( — Ij 2 ^ Ul - Un 



(7) 
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for an even n and 

I -I \ ( o I 

^ (2n - l)(2n - 3) ■ ■ ■ + 2) ■ ■ 

for an odd n. In T^j,^ ■.•!/„ is defined as 

where g'^ = g^'^ — v^v" and = c?^ + ^fi'^^- It is not difficult to check (H) is orthogonal 
to other states in @) and (H). For instance, (^) is orthogonal to the second state of (|) and 
the first state of (|), because / d'^pd^kTr{C'~^^(^^j_^...^^C} = 0. For the second state of we 
can write its matrix element of the current as 

S{pi---pn) 

E E ■ ■ ■ aT'" + Bgr^gr'gT' ■ ■ ■ + ■■■), 

5(l/l---I/n) S{pl---pn) 

(9) 

where the term proportional to A vanishes, the rest are also zero because 'j^'j'^u^'^'"''" vanishes 
after contraction of any two indices pi. 
Similarly we choose 

= r'i--'^"^--""/i(x)j,,...,„(x), (10) 

as the baryonic current of the second state of (|^), where 

YPl-PnVl-V„ _ ^ ^Pl ^Jj^l^M2 1^2 ^Pnl^n {^^) 

S{i---n) 

The currents for the rest states can be obtained by omitting the 75 inserted in the two light 
quarks. It does not change the quantum number (except isospin) of the currents to insert p 
and Ip into q'^q, or replace insertion of 7^ by V^, or let Dj_ act on q^ . But these insertions 
probably affects the isospin character. For instance, in the case of n = 1, currents (|^) and 
(p!0|) only have 1 = state. After the insertion of p, they only have 1 = 1 state. We will go 
back to discuss this later. Giving a definite quantum number, we only consider the currents 
with the lowest dimension, because higher dimension currents make QCDSR more unstable. 
So we will not discuss the currents with insertion of Ip later. 

The baryonic currents also can be obtained by using the wave functions proposed in 0. 
The result is the same as above. 

Now we are going on the renormalization of these baryonic currents, which must be used 
in the QCD sum rule beyond correction as well as matching the matrix elements of 
these currents from HQET to the full QCD. The one-loop renormalization of these baryonic 
currents without the covariant derivative has been considered in |^. When we consider 
the renormalization of the currents with one or more derivative D^, the currents with the 
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same quantum number and dimension will mix with each other. For convenience we firstly 
consider a symmetric tensor current, which is transverse, traceless and only has one covariant 
derivative D^, i.e. 

jT = grCT(Dl7l + /^l7l - \9Tmq2h. (12) 

Obviously the renormalization of ji^ is the same as J± in and (0) in the heavy quark 
limit. The complete set of operators which may mix with (0) is 

jT = q(CT{DlY^ + Dir^-lg^^ip)q,h 

jT = q(CT0W^ + 'Dl^^^-lgr'^^)q,h (13) 

where T could be 1 75, p or 75/; and q^^j_ = (D'^q)^ . j^^ is trivial because it differs from 
ii^ + by a total derivative. The nuisance operators such as 

qlCT{A^^Y^ + Al^'i-^-g^^4)q,h 

do not appear because they do not vanish by the equations of the motion. 
Then the bare current [ii"] can be expressed as 

[jT] = Z,Zl (Zijr + Z^jT + ^3jT), (14) 

where Zq = 1 — Cp-f^^ and Z/^ = 1 + 2Ci?^ are the renormalization coefficients of the light 
quark and the heavy quark resp. and = 1,3) are the renormalized currents. In order to 
determine the renormalization coefficients Zj(i = 1, 3), we need to insert the current into IPI 
diagrams and extract the ultraviolet divergence. We find that the divergence associated with 
the three-quark vertex are sufficient to determine all counterterms. The relevant Feynman 
diagrams are shown in Fig.l. In the Feynman gauge, we obtain 

Zi = l + [(f + 3/2)Cb-C7^]|^ 

Z2 = i^ + l)CBf^, (15) 
Z3 = 0, 

where -ff,TY±7^ = HTY±, Cb = ^,Cf = and e is defined as D = 4 - 2e in the 

MS-scheme. The renormalization of ^2"^ is similar as that of jj"^. Therefore ji'^ and j!^'^ can 
construct a renormalization-invariant space. The one-loop renormalization invariant currents 
can be chosen as 

J± — Ji ^ 32 ■ 

They satisfy the equation 

[jr] = ^,^l(^i±^2)jr. (16) 

The current j'^ corresponds to the state / = 1 for F = 1, or 75/; or / = for F = 75. j'^ 
corresponds to / = state for F = 1, or 75/; or / = 1 for F = 75. 
The anomalous dimensions are obtained by 

dlog{Zi ± Z2) 



7± 



dlog^ 



27r' 
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where /i is the renormahzation point. The relation between HQETs currents and the corre- 
sponding ones in the full QCD in the leading logarithmic approximation is 

Jqcd = C{ii)Jhqet{.Ij)^ (18) 

where the matching condition is 

r9(M) , i±(.a) 

C±iii) = C±(Mft)e^«('^^'^' ^ ^^'^ (19) 



and where C±{Mh) = 1 + 0{as{Mh)). It is more complicated to get the full matching 
condition. The relevant references could be seen in We do not discuss this case. 

(P^D can be extended to a general case. For the rank n tensor current with n — 1 covariant 
derivative D , the complete set of operators is 
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( — 1)^ V— \ * Ul---Un 

+ (2n-l)(2n-3)---(n + l) /^'^^^^ ' " " ^^'^"---^"--f ^)^(^ = 1'^) 

(20) 

for an even n and 

_ (rpi 1 rpi 1/ -L . . . 

Jvx---Un \-''Vi--Un 2n — 1 9 -^l^l'^2 UU:y ■ -Un ' ' ' ' 

+ (2n-l)(2n-3)---(n + 2) ^^^^^ ^^^^^^^^ ' " " = 1'^) 

(21) 

for an odd n, where 

= 9i(^rCT ^±.i---^±..-.^±.r--^±.„-i7±.„g2(x). (22) 

S(y\---Vn) 

The bare current [jf^ could be written as 

[jf -^1 = ^rf"'"- (23) 

fc=l,n 

After calculating the Feynman diagrams similar to Fig.l, we obtain 

1 1 _ m + 1 ^2 ^j^^ 



~' k\{n-l + l)\ 2''^A^e 

(24) 

for i = k, 

7 - f V^r- (^-OK^-/-l)!(^-fc + l)!(fc-/-l)! g^^ 

i(fc-^-2)!(n-fc + l)! a. ^ ^ 

^^'^-^ (n-. + l)! J 4^^ 
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for i < k and 



7 _ ^'x^^rm r^-m-iri {m - - I - l)\{n - k + l)\{k - I - 1)\ ^ 

^of^o t\{n-l + l)\ 2 (26) 
, ^n~i {-i-k-iy.kl as 
i\ J 4^ 

for i > k. The renormalization-invariant currents can be obtained through the diagonization 
of this matrix. We do not find a simple way to diagonize it. However, it should not be 
difficult to get the numerical solution for a special n. 



2 Two-point correlator 

The basic idea of the QCD sum rule is to consider the two-point correlator 

n^....M„,..-.„(fc) = / rf'xe^'=-^^(0|j;^...^Jx)J;^^...,jO) (27) 

for a region of k in which one can incorporate the asymptotic freedom property of QCD via 
the operator product expansion (OPE), and then relate it to the hardonic matrix elements 
via the dispersion relation. The currents ^^^...^^(a;) and Jl,^...^^{x) were defined in (^ or 
(p!0|). They have the same quantum number and the same F defined in (^). We do not 
discuss the case that these two currents have different F, because the perturbation part of 
such correlators is zero. The extra indices i,j were defined in (p2|) . 

In order to analyze the construction of the correlator ( P7D up to the order 1/Mh, we write 
down the effect Lagrangian of the heavy quark up to the order 1/Mh 

Leff = Kiv ■ Dh^ + -^{h^iyiDYK - \hvgG ^^a^^ . (28) 
ZlVlfi A 

In the leading order of HQET, the the spin of the heavy quark decouples with the light 
freedom, therefore (|27|) can be expressed as 



iio \^) - iiolt^j-'ai-Q„,/3i-/3n-'- ± J- ± ) l^yj 

where /ai...a„,/3i.../3„ is transverse, symmetric and traceless for indices a\ - ■ - an and Pi - ■ ■ Pn 

resp.. Since is proportional to via the heavy quark propagator J dtS{x — vt), only 
one scalar argument u = k ■ v survives in (pOl) . Joi...a„,/3i.../3„ is only composed of and g^u 
and must be written as 

2 

5(ai-a„)S{/3i-/3„) ^ 

-2r 



+ ■ ■ ■ + (2n - l)(2n - 3) ■ ■ ■ (n + if^"^"^9±P,P2 " ■ ■ 9±a,,.,a„9±p„.,pJ 



(30) 
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for an even n and 

2 

Iai—a„,l3i---f3 

n 

5(ai-a„,)5(/3i-/3„,) 

(-2)— 

"I" ■ ■ ■ "I" ^2j2 — l)(2n — 3) ■ • • (n + 2)^"*""^"^^^^^^^ ' ' ' 9-i-an-2a„-i9±Pn-2Pn-i9-La„P„) 

(31) 

for an odd n. 

To the order of l/Mh, the term j^hjj{iD)^hv does not break the spin symmetry but 
could produce a new scalar argument fc^ in the correlator p7|). For convenience, we choose 
k paralleling to v. Then we also can write its contribution as 

la l^^J = ^^iilaltUJ-'Qi-a„,/3l-/3n-'- ± ~2~ ' ^ ' 

The chromo-magnetic term -^^^KgC^^aj^hy breaks spin symmetry and causes the mass- 
split of the heavy baryon doublets. We can write its contribution as 

T-r/il---Ai„,l/l---I^„/ \ _ -T-r / \p/il---Mn,Ol--an^_J^ ^ + /' -fSvl---Vn,/3l---/3n ( nO\ 

where -Fjj,Q-i...a„,/3i.../3„ is transverse, symmetric and traceless for indices ai ■ ■ ■ a„ and Pi - ■ ■ Pn 
resp.. Since -Fjj,ai---a„,/3i - /3„ also only depends on and (7^,^, we can write down 



5(ai---an) 5(/3i-/3„) 

2 

5'±/3i/32 5'±«3/33 ■ ■ ■ 5'±a„-i/3„-i H i- 



2n - 1 

n-l 
2 



-2)^ 



(2n-l)(2n-3)---(n + 2) 

(34) 



for an odd n and 



-2/3n-2 

2 

~2n - l^-L"i"2fi'-L/3i/325'±a3/?3 ■ ■ ■ fl'±a„-2/3n-2 ^ 

(-2)— . 
(2n - l)(2n - 3) ■ ■ ■ (n + 3)^^«i«2^^/3i/32 " ■ ■ ^±a„_3«„-2^±/3„-3/3„-2 J 

(35) 



for an even n. 

It is not difficult to check that 



2 ' 



J-+ —<^±^ij,ai--an,Pi--l3n — -^nlai...a„^i3^...i3^L , 

pMi-Mn,ai-anlj^ iir. (36) 



i(n + i)/„,...„„,^,...^„r 



Mr 
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This is consistent with the naive expectation based on the quantum mechanics. 

With these formulas in our hand, we begin to construct the QCD sum rule for masses 
of the heavy baryons. At first, the imaginary part (^) can be obtained by insertion of the 
physical states. On the assumption that only one resonance and continuum contribute to 
the correlator (|27D , we can write 



J (T'xe''^''' {Q\J^^...^,^{x)\B^j^i) {B^j^i\Jy^...yJ\Q) + continnum 



-continuum. 



where M^_^i is the mass of the heavy baryon, A is defined as 



(37) 



M. + A + 0(_ 



and 



In (|38|) M is a I spinor and eai-a„ is a symmetric polarization tensor and satisfies 



(38) 



(39) 



Tj}'"'^"'"^'"'^" eai-a„u can serve as a n ± | spinor]^. There are 2n + 1 independent ea^...^^ 
which are normalized as 

2n+l 

E i * 



-/3l-/3n 



(40) 



i=l 



Then the sum rule can be established. The correlator (j^7\} obey the dispersion relation 



n±(a;) = - 



1 f ImU±{uj')duj' f p±{uj')duo' 



TT 



iu' — LU 



where the spectral density p±{uj) = ^Imll±{u) can be obtained from (|29|) , 
(p6|) (via OPE in a proper region of u) on the one hand. 



n±(cu) = no(cu) + 2ig-(n,„(cu) + a±Uib{u;)), 



a+ = n. 



fn + 1). 



(41) 
3D and 

(42) 



On the other hand, 



p^(uj') = tt6{uj - A)|/„^i pM^^t^ + continuum 



is obtained from (37). Thus, the baryonic mass as well as its matrix element can be extracted 



from (^ll). One may note that there are two sources of uncertainty in (0). At first, we have 
to truncate the series of OPE. Thus, the contributions from higher dimension operators raise 
the main uncertainty of QCDSR. Another uncertainty is from the contributions of heavier 
resonances and continuum. In order to suppress these uncertainties, the Borel transformation 



is applied on the both side of (41). Then we obtain 

'rf^e-^p±(^) = /„^./:^.M-ye 



(43) 



where M is the Borel parameter and the upper bound of the integral s is used to remove the 
rest contributions from the continuum. The mass parameter A can be obtained by taking 
the partial derivative respect to 1/M on both sides of P3|) 



Jo dujp±[uj)e M 

The Borel parameter should be chosen in a region that K± is not sensitive to it. In graphic 
language, there should be a plateau in the plot of M-A±. 



3 Sum rule for (3/2,5/2) 



For a special case, let us consider the heavy baryon doublet (3/2,5/2). As a preliminary 
calculation, the radiative corrections are not taken in account. So the heavy quark propagator 
to the order 1/M/j in the coordinate space can be expressed as 



dt l + rr^(-d^ -2qG,„.x^'d'' + ^(a.G^) -^^^a^"^^) 5ix-vt), 

(45) 

where we have used the fixed point gauge x'^yl^(x) = 0. The two-point correlators then are 



where 



x] 



-I [ TT{T/\^,^,{x)Si{x)/\p,pM^Sf{x)} 

+ (_l)/+lTr{5f(x)K,a,(x)^^5Kx)^^,;3,(0)f}]^^^^^'"^"^5,(x)n^'^^'^^^^ 

(46) 



A 



T = 7or^7o, = C^^r-^C, Sf{x) = C~^Sf{x)C, I is the isospin of the current and the light 
propagator Si{x) in external gluonic field was given in [Q. The correlators in momentum 
space related to (E6|) via Fourier Transformation 



n 



X] 



(47) 



Because we set paralleling to f^, by using the technique of the partial integral, we can 
write d^ ) in the form 



/ d'^xe'^^m^ 



(48) 



We analytically continue correlators from t > to imaginary t = —it. Then the spectral 
density p±{oj) are related to n±(r) by the Laplace transformation [Q 



ZTT J a—ioD 



(49) 
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Fig. 2 shows the Feynman diagrams that we calculate. The corresponding results are read 

as 

= 630^ - 864^^^"^^^' 

for r = 1 with / = and F = 75 with 1=1, 

= 630^ - 432^^^^'^' 

^ 140|r4 4327r3^ 

Pib{'^) = -——:{asG'^) 



for F = /> and F = 75/; with 1=0, 



p^^(^) = ^.^(a^G^), (52) 



for F = 75 with 1=0 and F = 1 with / = 1, 
for F =/>, 75/; with 1=1. 

It is interesting to compare (^) with (|5l|) . ( ^OD shows that / = state of the current 
with r = p and F = T^jb are degenerate up to the order we consider. However, from (|5T[), 
/ = state of the current with F = 1 is degenerate with 1 = 1 state of the current with 
F = F5 (not naively expected 1=0 state). It is another example that the insertion of p affects 
isospin. The degeneration of (^)-(|53|) will disappear when the light quark mass or the quark 
condensate is taken into account, i.e., the chiral symmetry is broken. Since our calculation 
in this paper is very preliminary, we leave the calculation of the quark condensates in our 
next paper. 

Now let us proceed to the determination of the heavy baryon masses. The mass of 
the heavy baryon can be obtained by (|4^) . In order to distinguish the mass split for the 
doublet(3/2,5/2), we expand the masses as 
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The sum rules for Aq and A^ are shown in Fig. 3 for the case (pOD, where we use the 
threshold s = 2GeV and the gluonic condensate {agC^) = 0.07(jeV"'^|TO|| [plj] . Fig. 4 shows 
the dependence of the Aq on s, where it is shown that s < Aq when s = 1.7GeV. The Aq 
is not very sensitive to s around s = l.SGeV. It goes down when (asG^) becomes smaller. 
The sum rule for (|5l| ) is very similar to the case (^). It gives the same mass prediction 
within the uncertainty of these two sum rules. The sum rule of (^2|) is much different from 
that of (0) and (0). Fig. 5 shows that the kinetic term is negative when s > l.GGeV. 
We do not find find a stable region for the threshold s. Fig.6 shows the sum rule of (|52D at 
s = lAGeV. The situation of is similar to f^2t). Thus we obtain 



(55) 



Ao = 1.8 ± 0.2GeV, Ak = 0.5 ±0.5GeV^, Aa/ = 0.05 ± 0.05Ge1/l 
for the sum rules ([50|) and (^T]) and 



Ao = 1.2 ± 0.2GeV", A^ = 0.06 ± OmGeV^, Am = 0.03 ± 0.005GeV^ (56) 

for the sum rules ([52|) and (|53D . The error bars are from varying the parameters in a proper 
region. Am is very sensitive to {asG"^). When {asG^) = O.OiGeV, it becomes unbelievably 



small. This is similar to the case in [Q. The reason may be that the mass split mainly is 
from the internal gluon exchange which we did not consider here. Ax/^GeV is much smaller 
that the masses of b and c quarks. Therefore the expansion of l/M^ stands well. Using the 
quark masses |l^ 

rrih = A.SGeV, m, = lAGeV, (57) 



we predict the masses of the baryon doublet (3/2,5/2) are 

Mb = 6.6 ± 0.2GeV, = 3.4 ± OAGeV. 



for (|5|) and (P) and 



Mh = 6.0 ± 0.2GeV, 



Mc = 2.8 ± 0.3^6^^, 



(5J 



(59) 



for (H) and (H). 

Finally, let us give a brief summary. We have discussed the baryonic currents with 
arbitrary quantum numbers as well as their one-loop renormalization. Using the obtained 
currents, we have analyzed their two-point correlators. For a special case, we did the QCD 
sum rule for the masses of the heavy baryon doublet (3/2,5/2) up to the order l/M^. Because 
we did not take account of the light quark mass and the quark condensates, we could not 
distinguish the two-point correlators of the currents with different parity. Also because we 
did not take account of the internal gluon exchange, the mass split of the heavy baryon 
doublet is too small. We will include these discussions in our next paper. 
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Figure captions 



Fig.l Feynman diagrams for the renormalization of the baryonic currents. 

Fig.2 Feynman diagrams for the two-point correlator. The smaller dots stand for the 
current vertices, the bigger dots stand for the l/M/^-order interaction vertices. 

Fig.3 The Sum rule for (|50D. Dashed line gives Aq versus Borel variable M, solid line is 
Kk/^GcV versus Borel variable M and dotted line is Ajvf/lGeV^ versus Borel variable M. 

Fig.4 The dependence of Aq on s. Dashed line gives the result for s = l.SGeV, solid line 
is for s = 1.7GeV and dotted line is for s = 2GeV. 

Fig.5 The dependence of Ak/IGcV on s for (^21) . Dashed line gives the result for s = 
lAGeV, solid line is for s = 1.2GeV and dotted line is for s = l.GGeV. 

Fig.6 The Sum rule for (|52D. Dashed line gives Aq versus Borel variable M, solid line is 
Ax/lGeV versus Borel variable M and dotted line Ajv//lGeV^ versus Borel variable M. 
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